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On the example of the spherical model we study, as a function of the temperature T, the 
behavior of the Casimir force in 0(n) systems with a diffuse interface and slab geometry 
oo"- 1 x L, where 2 < d < 4 is the dimensionality of the system. We consider a system 
with nearest-neighbor anisotropic interaction constants Jm parallel to the film and J± across 
it. The model represents the n — > oo limit of 0(n) models with antiperiodic boundary con- 
ditions applied across the finite dimension L of the film. We observe that the Casimir am- 
plitude Acasimir(d| Jx% J\\) of the anisotropic d- dimensional system is related to that one of the 

isotropic system Acasimir(d) via Ac aS imir(d| J±, J\\) = (J±/J\\) ( ' d ~ 1 ^ 2 Ac as imir(«0- For d = 3 we 
find the exact Casimir amplitude Acasimir = [CI2 (t/3) /3 — £(3)/(67r)] (J±/J\\), as well as the ex- 
act scaling functions of the Casimir force and of the helicity modulus Y(T, L). We obtain that 
/3 C T(T C , L) = (2/tt 2 ) [Cl 2 (vr/3) /3 + 7C(3)/(30tt)] (J±/J\\) where T c is the critical temperature 
of the bulk system. We find that the effect of the helicity is thus strong that the Casimir force is 
repulsive in the whole temperature region. 



PACS numbers: 05.20.-y, 05.50.+q, 75.10.Hk 



I. INTRODUCTION 



The excess free energy due to the finite-size contribu- 
tions to the free energy of a system with a film geometry 
characterizes a fluctuation-mediated interaction which is 
termed the Casimir force, or, in the case of a fluid con- 
fined between two parallel walls - also the solvation force 
(or the disjoining pressure). The force is named so af- 
ter the Dutch physicist Hendrik B. G. Casimir who in 
1948 [1] first noticed that when two metallic perfectly 
conducting uncharged plates face each other in vacuum 
at zero temperature the restriction and the modification 
of the zero-point vacuum fluctuations of the electromag- 
netic field between the two parallel plates lead to a de- 
pendence of the energy of the system on the distance L 
between the plates and, thus, to a force between them 
which turns out to be attractive. The above is the so- 
called classical (actually quantum mechanical) Casimir 
effect. When the fluctuating medium is not a vacuum, 
but a thermodynamic system, say fluid, near its bulk crit- 
ical point T c one arrives at the so-called thermodynamic 
Casimir effect that has been predicted by M. E. Fisher 
and P. G. de Gennes [2] in 1978 and which has been a 
subject of intensive theoretical and experimental studies 
afterwards [3-25]. 

We remind that for an 0(n),n > 1 model of a d- 
dimensional system with a temperature T and geometry 
oo d_1 x L the Casimir force is defined by [3], [5] 



where f^(T,L) is the excess free energy 

f£>(T,L)=f( T \T,L)-Lf b (T), 



(1.2) 



and the superscript r denotes the dependence on the 
boundary conditions. Here f^(T,L) is the full free en- 
ergy per unit area of such a system under boundary con- 
ditions t and fb is the bulk free energy density. It is 
believed that if the boundary conditions r are the same 
at the both bounding the system surfaces F { 



(r) 

Casimir 



F, 



(t) 



Casimir 



(T,L) 



df£\T,L) 
dL 



(1.1) 



will 

be negative. In the case of a fluid confined between iden- 
tical walls this implies that then the net force between 
the plates will be attractive for large separations. If the 
boundary conditions are essentially different at the both 
confining the system surface planes (e.g. one of the sur- 
faces prefer the liquid phase of the fluid while the other- 
prefers the gas phase) the Casimir force is expected to 
be positive in the whole region of the thermodynamic 
parameters, i.e. then the net force between the plates 
will be repulsive. 

In the current article we will investigate the behavior 
of the Casimir fore in systems with diffuse interface. As 
a realization of such systems one can think of about the 
reaction of 0(n) models with n > 2 to some helical exter- 
nal field which reaction can be characterized in terms of 
some helicity modulus T or, in case of a magnetic mate- 
rials, of Bloch walls between the domains of the magnet. 
Heuristically, the helicity modulus is the analog of the 
interface tension for 0(n)-symmetric systems. The sim- 
plest theoretical model of a system with a diffuse interface 
is the 0(n) model with antiperiodic, i.e. r = a, bound- 
ary conditions and short-ranged interactions. According 
to the standard finite-size scaling theory (see, e.g., [5, 26] 
for a general review) one expects that near the critical 
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temperature T c (of the corresponding bulk, i.e. L = oo 

(a) 

Casimir 



system) the behavior of ^casimir wm ^ e & Yen by 



(1.3) 



^:Li r (r,i)=i- d 4 a a ) sim i r (^), 

while that one of the full free energy is 

P.f^{T,L) = L-^xf\x t ), (1.4) 

where x t = a t tL x l v is the temperature scaling variable 
with t = (T — T c ) /T c being the reduced temperature, 
a t is a nonuniversal scaling factor, while X^ simil and 
xj a ' ) are universal (geometry dependent) scaling func- 
tion and v is corresponding (universal) scaling exponent 
that characterizes the temperature divergence of the cor- 
relation length when one approaches the bulk critical 
temperature from above, i.e. £(,(i — > + ) ~ £o~i _ ". The 

scaling functions v ^ 
relation 



Casimir an d ^/ are related via the 



4 a), 



4t mir (**) = (d-l)^^)-^* ^^r^)- (1-5) 

The value of at the critical point is known as the 

Casimir amplitude A< a ), i.e. A< a ) = X { f a) (x t = 0). On 
its turn, the excess free energy under antiperiodic condi- 
tions f$ can be related to the one of the same system 
under periodic boundary conditions f$ via the finite- 
size helicity modulus T(T, L) [27] 

/&>(T,L) = f£(T,L) + ^T(T,L), (1.6) 

where T(T) = linii,^ T(T, L) with T(T) > 0. For 
the behavior of T(T, L) near T c the standard finite-size 
scaling theory states that 



f3T(T,L)=L-^X r (x t ), 



(1.7) 



where Xx is universal scaling function. Actually, when 
d = 3, a modification of Eq. (1.7) has been suggested in 
[28] by Privman, who supposed the possibility of appear- 
ance of "resonant" logarithmic term due to the mutual 
influence of the regular and singular contributions in the 
helicity modulus 

/3T(T,L) = L- 1 \Xr(xt)+ujHL/a)\+^(T)L- 1 + --- , 

(1.8) 

where u> is an universal amplitude, while <1>(T) is a regular 
at T c function and a is some characteristic microscopic 
length scales (e.g., the distance between the molecules 
of the correlated fluid, or the lattice spacing). The va- 
lidity of this hypothesis has been checked in [27] on the 
example of the exactly solvable mean-spherical model. 
No logarithmic corrections of the type predicted in (1.8) 
have been found. Let us recall that in the case of super- 
fluids (n = 2, d = 3) the helicity modulus T is propor- 
tional [29] to the superfluid density fraction g, namely 



q = (m/h) 2 T(T) with m being the mass of the helium 
atom, and is directly measurable (for experiments mea- 
suring q in thin films of 4 He see, e.g., Rcfs. [30] and 
[31]). In fact, (1.8) was proposed in [28] as an attempt 
to improve the fit of the experimental data. It turns 
out, however, that the overall fit of the data is improved 
only in a very limited way, provided one insists on the 
bulk value of v in the scaling variable x t . The scaling 
" data collapse" technique works well if one takes v as an 
adjustable parameter not necessarily equal to the correla- 
tion length exponent. It also should be emphasized that 
one could expect additional complexity in the behavior 
of the finite-size scaling function of the helicity modulus 
in the case of superfluid transitions in a film geometry; 
nevertheless, the analysis of the experimental data shows 
no clear singularities or a jump in the finite-size scaling 
function [30, 31]. 

According to all the accumulated analytical and nu- 
merical evidences, see e.g. [5], [26] and references cited 
therein, when a: t > 1 both the excess free energy and 
the Casimir force under both periodic and antiperiodic 
boundary conditions in systems with short-ranged inter- 
actions is expected to tend to zero in an exponcntial-in-L 
way. This is consistent with T(T) = for T >T C . When 
x t — > — oo the same quantities tend to zero in an power- 
law-in-L way. This slow algebraic decay of / ex (and of 
-Fbasimir) is, of course, associated with the existence of 
soft modes in the system (spin waves) when T <T C and 
in the absence of an ordering external field destroying 
the 0(n) symmetry. This, in turn, will lead to a much 
greater (in comparison with the Ising-like case) Casimir 
(solvation) force when T < T c in 0(n) models. With 
respect to the Casimir force the last has not only been 
predicted theoretically, but has been also observed exper- 
imentally [10, 14] and, relatively recently, confirmed in a 
model study of the XY model numerically via Monte 
Carlo simulations [22, 23]. The considered systems do 
not posses, however, a diffuse interface. When such am 
interface is present and T < T c from Eq. (1.6) it is easy 
to see that 
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Casimir 



(T<T C ,L> 1) ~ 1-tt 2 T(T)L- 2 . (1.9) 



Since T(T) > the last implies that the force will be 
repulsive and much stronger, of the order of L~ 2 , than 
in systems with no diffuse interface where it is either of 
the order of L~ d , or smaller. 

Since we consider film geometry, it is natural to al- 
low for an anisotropy of the interactions in the system 
which reflects this geometry. To that aim we will take 
the interaction constant in the Hamiltonian along the 
surface, say J\\ , to be different from the one perpendic- 
ular to the film, say J±. Since such anisotropy does not 
change the universality class of the bulk system one might 
naively expect that the scaling functions of the finite sys- 
tem X^ simir , xj a ' and Xf will be the same as for the 
isotropic system. Recently it has been argued, however, 
see Refs. [32, 33], that this is not true and that one shall 
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expect these functions to be nonuniversal and depending 
on the ratio J±/Ju. It has been shown [32, 33] that the 
main reason for this state of affairs is the need of a gen- 
eralization of the standard hyperuniversality hypothesis 
[34-39]. According to it, if fb, s ing(T) is the singular part 
of the bulk free energy density /& normalized per ksT 
and £(T) is the bulk two-point correlation length in the 
isotropic system, then 

lim f b , sing (T)[£(T)] d = Q, (1.10) 

where Q is a universal constant that characterizes the 
corresponding universality class. If now fb,sing(T\Jj_, Jy) 
is the corresponding free energy in the anisotropic film 
system with £ii being the correlation length along the 
system surface and £j_ the one perpendicular to it, then 
the generalized hyperuniversality hypothesis states that 

lim J bMng (T\J ± , J ]] )[^(T)} d - 1 a(T) = Q, (l.H) 

T->T+ 

with Q being the same universal quantity as in the 
isotropic case. Note that the new hypothesis involves 
two different correlation lengths, characterized by two 
different correlation length amplitudes, while the stan- 
dard hypothesis deals with only one correlation length. 
Note also that the validity of (1.10) is one of the main 
prerequisites for arguing the validity of the scaling hy- 
pothesis (1.4) by Privman and Fisher [40]. It is, however, 
possible to relate the scaling functions of the anisotropic 
to that one of the isotropic system. Indeed, choosing the 
isotropic system to be such that its correlation length is 
equal to, say, £j_ from (1.10) and (1.11) one obtains that 



/fc,sing(r| Jj_, J||) 



and, thus one arrives at 



fb,sing(T), T > , 

(1.12) 



Xf\x t \j^^) 



X<f\x t ), 



(1.13) 



where £y and £j_ are the correlation lengths in the 

anisotropic system while X^\xt) is the universal scal- 
ing function of the isotropic one. Of course, (1.12) and 
(1.13) shall be considered only as plausible hypotheses 
which validity has to be verified. Note that, if valid, Eq. 
(1.13) implies a relation of the Casimir amplitudes in the 
anisotropic and isotropic system 

^Casimir Casimir (d). (1.14) 

In the current article on the example of the exactly 
solvable mean spherical model with 2 < d < 4 we will 
demonstrate that in the anisotropic system with a diffuse 
interface the scaling function X^ simil , xj a ' and Xr in- 
deed depend, in addition on the scaling variable x t , also 
on the ratio Jj_/J\\. This will lead, e.g., to nonuniver- 
sality of the Casimir amplitudes in such systems which 



are, however, simply related to the ones of the isotropic 
system via the relation (1.14). We will determine the 
explicit form of the scaling function of the free energy, 
Casimir force and of the finite-size hclicity modulus. For 
the case d = 3 in the isotropic system we will find the 
universal values of these quantities at the critical point 
T c of the bulk system. We will also consider the case 
when the nearest neighbor interaction Jy along the film 
might be different from the one in orthogonal direction 
J±. 

The structure of the article is as follows. In Section 
II we define the model under consideration and provide 
some basic expressions needed for its treatment. The 
results for the finite size behavior of the free energy and 
of the Casimir force are presented in Section III, where 
in subsection III A we present our general results for 2 < 
d < 4, while in subsection IIIB the explicit results for 
the important case of d = 3 are given. Our findings 
about the behavior of the helicity modulus are contained 
in Section IV. The article closes with a discussion and 
concluding remarks given in Section V. Some technical 
details and results needed in the main text are derived 
in Appendixes A and B. 



II. THE SPHERICAL MODEL 

As stated above, we will study the finite-size behavior 
of an anisotropic system with a diffuse interface on the ex- 
ample of a spherical model embedded on a (i-dimensional 
hypercubic lattice £gZ d , where C = L\ x Li x • • • Ld- 
Let Li = Nidi,i = 1, • • • ,d, where Ni is the number of 
spins and a* is the lattice constant along the axis i with 



&i being a unit vector along that axis, i.e. Bi.ej 



Si, ■ 



With each lattice site r one associates a real- valued spin 
variable S r which obeys the constraint 



rEC 



(2.1) 



where N = NiN 2 ■ ■ ■ is the total number of spins in 
the system. The average in (2.1) is with respect to the 
Hamiltonian of the model 



-i/35> r J(r,r')Sw+ S ]TS r 2 . 



(2.2) 



In the current article we will consider only the case of 
nearest-neighbor interactions, i.e. we take J(r, r') = 
J(\r — r'\) = Ji, if r — r' = iojej, i = 1, • • • , d, and 
J(r,r') = otherwise. Explicitly, one has J(r,r') = 

Z)»=i Md( r - r' - aiei) + S(r - r' + ajej)]. Let peri- 
odic boundary conditions are applied across directions 
e,, i = 1, • • • (d — 1), while antiperiodic boundary condi- 
tions, responsible for the creation of a diffuse interface 
within the system, are applied across e^. Generalizing 
for the considered here anisotropic case the results of 
[5, 7, 27, 41, 42] pertinent to an isotropic model, it can be 
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shown that the free energy of the model (per unit spin) 
is given by [43] 



+ sup £ \n\ S -Uj^(k) 



(2.3) 



feeBZi 



where N = (N u N 2 , ■ ■ ■ N d ), J = (J u J 2 , ■ ■ ■ , J d ), 
j( a ) (fe) is the Fourier transform of the interaction J, i.e. 



j( a )(fc) =J2 J{r)e lk r , 



(2.4) 



Jmax = max fc j( a )(fc), and the wave vector k = 
{hi, k 2 , ■ ■ ■ , k d } <G BZi is with components fc, = 2imi/Li, 
where rij = 0, • • • , Ni — 1, i = 1, • • • , (d — 1), while 
fed = 27r(n d + l/2)/L d with n d = 0, • • • , N d - 1. Thus, 
explicitly one has 



d-l 



j»(fc) =2^J, 



27TO^ 

"/vT 



2 Jd cos 



7r(2n d + 1) 



(a) 



(2.5) 



2 £?=i -A + 2J d cos(7r/iV d ) = Jo - 2J d [l 



and J, 

cos(7r/iVd)], with J = 2^" =1 Jj. Note that the ground 
state energy Jmax depends on N d and is twofold degener- 
ate - it is reached for both k\ = k 2 = ■ ■ ■ = k d -i = k d = 
and fci = k 2 = ■ ■ ■ = k d -\ = 0, k d = N d — 1. The equation 
(2.1) for the spherical field s reads 

^ E T77T77T77T = 1- (2-6) 



fceBZ 



!-ij8»(fc) 



We will be mainly interested in determination of the 
Casimir force and the helicity modulus within the con- 
sidered model in a film geometry. For that aim let us 
take Ji = J 2 = • • • = J d -i = Jy , N x = N 2 = ■ ■ ■ = 
N d -\ = iV|| , Jd = Jl, N d = N± and to perform the limit 
7V|| — > oo, i.e. to consider a system with a film geometry 
in which all the interactions in directions parallel to the 
film surface are equal (to Jy) but possible different from 
the interaction in the direction perpendicular to the sur- 
face (which is Jj_). Then, Eqs. (2.3), and (2.6) become 



f(«) 

"max 



PfW(l3,N ± \d,J) = llng-ix j 



(2.7) 



■ sup \u( a \w,N ± \d, J) - ^Kw\, 
w>o 2 J 



N± kd A||£BZi 



K 



w + uj^){k h k d \d, J)' 



(2.8) 



correspondingly, where k = (k\\,k d ) with fe|| = 
(fci,fc 2 ,--- ,fcd-i), K = /3J , 

U ia) (w,N ± \d, J) 

1 _ /-(d-i) 



27V 



E 



In 



J- fed Jfc||eB Zl 



w; + w (a '(fc||, fc d |d, J) ,(2.9) 



with 



u>)(fc|d, J) - [Jit - J (a) (fe)l /Jo > 0, (2.10) 



,(d-l) d-l f2v „ 



(2.11) 



and we have replaced the spherical field s by another field 
w, defined as 



Here 



w = 2s/K-j£ljJ . (2.12) 



Jit = 2(d 1) J,, + 2J ± cos(tt/N ± ) (2.13) 



is the ground state energy of the finite system under an- 
tiperiodic boundary conditions, while 

Jo = 2{d- 1)J|| +2J ± (2.14) 

is the ground state energy of the infinite one and, thus, 

J^L/Jo = (d- 1)&|| + b ± cos(7r/iV ± ), (2.15) 

where 

d 

b± = J±/J2 Ji ( 2 - 16 ) 



i=l 



and 



6 n = J n/E J » 



(2.17) 



reflect the asymmetry in the interaction. 

Eqs. (2.7) - (2.12) provide the basis for the investiga- 
tion of the behavior of the Casimir force within mean- 
spherical model in the presence of a diffusive interface in 
the system. 



III. FINITE-SIZE BEHAVIOR OF THE FREE 
ENERGY AND THE CASIMIR FORCE 

A. General results for the case 2 < d < 4 

From Eq. (2.7) for the excess free energy 

0fg>(J3,N ± \d,b) = N ± \j3fW(J3,N ± \d,b)-l3f b (J3\d,b)] 
one obtains 



0fg>(J3,N ± \d, b) = N ± 



-\k (w - w b ) + (w, N ± \d, b) - U d (w b \b) 



where f b (/3\d,b) = lim iVj .- > oo /(/?, N ± \d,b), b = 
(fen , • • • , by, b), w = w(K, N±\d, b) is the solution of Eq. 



■5 



(2.8), and w b = w b (K\d,b) is the lmijv ± ^oo limit of 
w b (K\d,b), i.e. w b (K\d,b) = lim Wi ^oo w(K, N ± \d, b). 
As it is well known, see e.g. [5], for K < K c = Wd(0\b) 
the spherical filed w b is solution of the equation 



K = W d (w b \b), 



where, for w > 0, 

i r (d) 

W d {w\b) = - / 



(3.2) 



(3.3) 



w + uj(k\d, b) ' 

and w b = 0, when K > K c . In Eq. (3.1) U d (w\b) = 
linijv^^oo U^ a \w, Nj_\d, b) which limit, according to Eq. 
(2.9), reads 



1 

Ud{w\b) = - \n[w + w{k\d,b)\- 

z JfeeBZi 



(3.4) 



Note that it does not depend on the boundary conditions. 
Obviously, the only nontrivial N± dependence in fix 
stems from the size dependence of the spherical field w 
and from the asymptotic behavior of U^ a \w, Nj_\d, b) on 
N± for N± 3> 1. Let us now study these dependencies in 
detail. 

Using the identity 



In a = 



00 dx , 
— ( e ~ 



'o x 

one can rewrite Eq. (2.9) into the form 



U {a) (w,N ± \d, b) 



1 



dx 



e — e 



, rf— i 



(3.5) 



S^(xb ± )x 



(3.6) 



where 



JV-l 



n=0 



7r 7r(2n + 1 

— z I cos — — cos 



N 



N 



With the help of the identity 



(3.7) 



Sffi (z) = cxp z (l — cos 



7T 

N 



2S ( ^(z)-S ( *\z) 



where 



(3.8) 



S%\z 



) = ^E cx p[- z ( 1 - cos ^ 

n=0 1 V 



(3.9) 



the problem for determination of the asymptotic behav- 
ior of the sum S^\x) when N 1, which character- 
izes the antiperiodic boundary conditions, can be reduced 
to the determination of the asymptotic behavior of the 
sum sff (x) , which is pertinent to systems with periodic 
boundary conditions. It can be shown that [27] 



c+(a) 
N 



-(a) _ 



'N 



e- x I (x) - 



= exp 



z(! - cos— ) 



2 *<-> ( ™L 

TlX V X 



(3.10) 



x < N , 



where 



R^(x) = J2^ in(n+1 



)x 



n=l 



x/4 



n=l 



n=l 



v(x) = —= [l — erf (71- v 7 ^)] • 

V 47TX 



(3.11) 

(3.12) 
(3.13) 



In addition, with the help of the Poisson identity, one can 
easily check that the following equivalent representations 
of functions R( + \x) and R(~\x) are valid 

R( + \x) = le x e 2 (0,e- 4x )-l, (3.14) 



R(-)(x) = ^(-l)"e-" 2a; / 4 = i ]e A (0,e- x / 4 ) - 1 



n=l 



(3.15) 
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where 82 (x) and 64 (x) are the corresponding theta func- 
tions. 

If one insists on using only the second asymptote in 
Eq. (3.10) as the one valid for all x, see, e.g. Ref. [42], 
then the corresponding result for U^ a \w, Nj_\d, b) reads 



(w^N^b) = U d (w\b) 



r^x- 

Jo x 



( 47r )d/2 y ^ J 
d/2 e -y XR (-) 



(d-l)/2 



where 



w = w — b± 1 — cos — — 
N ± 



(3.16) 



(3.17) 



y = y - n , 

and, see Eq. (3.6), 



y=(2N 2 Jb A _)w (3.18) 



1 f°° cIt ( 

U d (w\b) = -J ™\e~ x - e-* x [e- xb -I (xb ± )] 



n d-l 



(3.19) 



Using the representation (3.19) it can be shown [17] that 
when w — > + one has 



U d (w\b) = U d (0\b) + \wW d (0\b) 



1 r(-d/2) 
"2(27r)^ntiV^ 



■w d ' 2 



(3.20) 



with the dots representing terms of higher order than 
those retained in the expression. FromEqs. (3.16), (3.20) 
and with the help of the representation (3.15), for the 
finite-size part U (w, N±\d,b) of the free energy in the 



limit w — > + and, thus y > - see Eq. (3.18), one 
obtains 



U [a) {w,N±\d, b) 



= U d (0\b) + h ± y W d (0\b)N^ 2 - N^ d (j± 



(d-l)/2 



2 (47r) d / 2 (27r) d / 2 



ra=l 



+0(iVI 4 ). 



(3.21) 



Then, from Eqs. (3.1), (3.20) and (3.21) for the excess 
free energy one derives the final result 

Pf^i^N^b) = N-^(^~Y I \x t {y-y b ) 

i r(-d/2) / d/2 d/2 x 

2 (4tt)<V 2 ^ ^ ) 

1y d ' 2 y, nB ^/ a (nVy) l 
W' 2 ^ J (nVy) d / 2 J' 

(3.22) 

where x t is the temperature dependents scaling variable 

(d-l)/2 



(Jfc-Jf) i/=l/(d-2), 

(3.23) 

? y fc = (2iY|/6 ± ) Wb , (3.24) 



with Wf, being the solution of the bulk spherical field 
equation (3.2). 

Let us now see what is the correct answer when the 
complete asymptotic behavior, as given in Eq. (3.10), is 
used for the determination of the excess free energy. 

Using the asymptotes given by Eq. (3.10) one obtains, 
see Appendix A, that: 



U^(w,N ± \d,b) = Ud (0\b) + \b ± (y-n 2 ) W d (0\b) - itfj* (V)^ ^ _2_|r(-rf/2) y d / 2 



+7r 2 r(i - rf/2)// 2 - 1 + 2V4^ [ C 

Jo 



f^--(d-l)/2.-i/a: 



1 + U 2 2 



1 



2V47TX 



(1 +7T 2 X) 



(3.25) 



where y = (2N]_/b±)w > 0. Then from Eqs. (3.1) and (3.20) for the excess free energy one obtains 



0fg>(l3,N ± \d,b) = N-V-V( b -± 



(d-l)/2 



-x t (y - 7T 2 - y b ) 



1 1 

"2 (47r) d / 2 



r (-d/2) (// 2 - yf 2 ) + 7r 2 r (1 - d/2) y^ 2 - 1 



+ I(y,d)\, (3.26) 
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where 



(4„-)(d-l)/2 



1 + tHx 



dxa; -(d+l)/2 e -^ 



1 



2V47TX 



(3.27) 



The expression (3.26) has to be compared with Eq. (3.22) 

fa) 

that follows when one uses as asymptote of S N , when 

N » 1, only the asymptote from Eq. (3.10) (see, 

e.g., Ref. [42]). As we see, (3.26) and (3.22) differ from 
each other. However, using the identity 



1 + RW 



(3.28) 



one can show that when d < 4 and y > it 2 



i(y,d) 



-y d/2 



_2_ n „%M) 



1 T (-d/2) 
"2 (47r) d / 2 



^/ 2 ) +7r ^//2-l 

(3.29) 



and, thus, expression (3.26) is equivalent to (3.22) for 
y > ir 2 . In the opposite case, when y < 7r 2 , one can 
use (3.26) or, cquivalently, the analytical continuation 
of (3.22). Therefore, although in the derivation of (3.22) 
the incomplete asymptotic behavior of sums involved has 
been used, which makes this derivation mathematically 
wrong, and the expansion (3.20) of the bulk quantities 
has been applied, which is valid only for y > n 2 , Eq. 
(3.22) is still valid and can be used for all y > since 
this equations is equivalent to (3.26) which is obtained 
when one follows the proper mathematical procedures. 

When \y — n 2 \ < A-k 2 one can provide a representation 
of the integral I(y, d) in terms of power series which is 
very convenient for analysis of its behavior for small val- 
ues of the argument y. The corresponding representation 
is derived in Appendix B, and reads 



IM = 2^7* [^-d/2)+yT(-d/2)] 

CO 

-nW^atftf-yr, (3.30) 



(d-2)/2 



m=0 



where the coefficients a™ are given by 



,(<*) 



{2 x - d - 2- 2m )T(m + 1=2) C(2to + 1 - d) 



7r 2m m | 

(3.31) 

From Eqs. (1.1), (3.22) and (3.26) for the Casimir force 



one obtains the two equivalent representations: 
/3F^(0,N ± \d,b) 



(d-l)/2 



-x t (y-y b ) - (d- 1) x 



2 (47^/2 V Vb ) + * 



(2n) d / 2 



(3.32) 



in the derivation of which we have used the identity 
d 



dy 



y^K^ay) = -ay^K^ay), 



(3.33) 



and 



f h \ (d-l)/2 f 1 

= N ^{t) \ rt (y-* 2 -y b )-(d-i)x 



l l 



| 2 (4^) d /2 
+ir 2 T(l-d/2)y d / 2 - 1 



(-d/2) (> 2 y d J 2 ) 



i(y,d) 



(3.34) 



In (3.32) and (3.34) the variables y (or y) and y b satisfy 
the spherical field equations (2.8) and (3.2), respectively. 
It can be easily shown that these two equations can be 
rewritten in a scaling form. In the geometry of a film and 
under antiperiodic boundary condition the equation for 
y reads 

1 r(l-d/2) d/2 _ ! - rf /2-! 2 

2 ' ~ (47T) d /2 V +V (27T) d /2 

^(-ir ^ , (3-35) 
^ (nVy) d /2 1 

which is equivalent to 

1 rq-d/2) d/2 _ 1 r(2-d/2) 

2 ' (47T) d / 2 ^ 2 d 7T d /2- 2 * 



while the corresponding equation for t/6 is 

i _ r(i-d/2) d/2.! 



( 47r )d/2 



2/6 



(3.36) 



(3.37) 



Eqs. (3.23), (3.32), (3.34), (3.35), and (3.37) demon- 
strate that the Casimir force in a system with anisotropic 
interaction can be written in the form 

PFtiimiM N ±\d, b) = N^ d (^j ( ]l X Casimil (x t ) 7 

(3.38) 
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where Acasimir is a universal scaling function, provided a 
suitable definition of the scaling variables, see Eq. (3.23), 
is used. Note that x t is of the form x t — a t (b) tL 1 ^ 
which means that all the effect of the anisotropy of 
the type considered can be incorporated in the factor 

(& ± /6||) (d ~ 1)/2 = (J ± /J||) (d ~ 1)/2 in front of the scaling 
function on the r.h.s. of Eq. (3.38) and in the nonuni- 
versal factor at that enters in the definition of the tem- 
perature scaling variable x t . Note that with respect to 
the Casimir amplitudes Eq. (3.38) leads to the following 
relation between the amplitudes in the anisotropic and 
isotropic systems 



^Casimir (d\J±,J\\) 



J\\J 



(d-l)/2 



A CaS imir(d). (3.39) 



Note also that, because of the universality, the value of 
the Casimir amplitude in the isotropic system does not 
depend on J = J± = Jm. In order to achieve a confor- 
mity with the relation (1.14) one needs only to determine 
the ratio in the anisotropic system. In fact, this 

already has been done in [17] with the result that 



Si 



(3.40) 



Inserting (3.40) into (3.39) one, indeed, immediately ob- 
tains (1.14). 



Results for the case d = 3 



asymptote 



D 




FIG. 1: The scaling function Xcasimir(a;t) of the Casimir force 
F climiAP, N ± \ d = 3, b) for d = 3. Note that X C asi m ir(:r t ) > 
for all x t . The asymptotic behavior of -Xcasimir(a^) for xt -C 
— 1 is given according to Eq. (3.51). 

Since d = 3 is of special importance we will present 
some explicit results for this case. With d — 3, the equa- 



tions (3.32) and (3.34) simplify to 

/^:L ir (/w<*=3,&) 



-3/2 3/2 

y ' - v b 



(3.41) 



and 



P4 a limiM^\d = 3,b) 



+ ^V2 + 27(2/)3) | 5 



3/2 3/2 

y ' -y b 



(3.42) 



while the equations (3.35) and (3.36) for y and y become 



and 

x t 



1 r~ 1 , 

27T 7T 



1 



— r- — 



1 + 7T X 



1 + e 



dx 



-Vv 



(3.43) 



2\/47ra; 



l + i? (+) (7T 2 x) 

(3.44) 



respectively. Eq. (3.43) can be explicitly solved in the 
form 



(3.45) 



\fy~ = 2 arccosh 
At T = T c , i.e. when Xt = 0, this solution simplifies to 

r- 2-K 

V^-±z— . (3.46) 

As it is well known [5] , the scaling form of the solution 
of Eq. (3.37) for yb for the infinite system with d = 3 is 



27TXi, x t > 

0, x t < 



(3.47) 



At T = T c with y b = 0, according to Eq. (3.47), and 
y from Eq. (3.46) one can from Eq. (3.41) obtain the 
Casimir amplitude in the form 



^Casimir 



J± 



-Im (Li 2 (^T)) 

3 V K >> 6tt 



, (3.48) 



which, using the relation Im(Li2(e l9 )) = Ci2(#) between 
the polylogarithm and the Clausen function (see, e.g., 
[44]) 



a 2 (0) = £ 



*;=i 



sin(fc^) 



(3.49) 
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can be written as 

asimir - 



J__ 

h 



(!) 



03) 

6ir 



- 0.274.14;-! ( 



(3.50) 



One can also determine the full temperature dependence 
of the Casimir force. For that aim, in Fig. 1 we present 
the scaling function Acasimir (%t ) of the Casimir force 



F, 



(a) 



Casimir 



(/?, iVj_ \d = 3, 6) as a function of the temperature 
scaling variable x t . We observe that Acasimir^t) > for 
all x t , i.e. the Casimir force under antiperiodic bound- 
ary conditions is always a repulsive force. Furthermore, 
from Eqs. (3.41), (3.43) and (3.47) it is easy to check 
that when xt S> 1, one has y,yt,^$> 1 which lead to the 
result that the scaling function ^Casimir(^t) decays expo- 
nentially fast to zero, while for xt <C — 1 one has y — > + , 
yb = and that 



Aca 



r(ar t ) 



-tt 2 x/4- C(3)/tt. 



(3.51) 



As we will see below, the last equation, together with 
Eqs. (3.38) and (4.7) - see below, lead to the conclusion 
that when T«T C the behavior of the Casimir force in 
systems with a diffuse interface in indeed given by Eq. 
(1.9). 



IV. HELICITY MODULUS 

A. General results for the case 2 < d < 4 

The concept of the helicity modulus was introduced by 
Fisher et al. [29]. Fundamentally, the helicity modulus 
is a measure of the response of the system to a helical 
or "phase-twisting" field. Alternatively, for an isotropic 
system with n-component order parameter, where n > 2, 
one can consider the helicity modulus to be the analogy 
of the surface tension or interfacial free energy between 
two phases in a system with a scalar, i.e. n = 1 order 
parameter (e.g., an Ising model). In other words - the 
helicity modulus is a measure of the increase of the energy 
of the system due to the existence of a diffuse interface 
within it. When in an 0(n), n > 2 such an interface 
is created by the application of antiperiodic boundary 
conditions the helicity modulus can be defined, e.g., as 
suggested in [27] 



T(p,N ± \d,b) 
2A ± 



/£> (/?, N ± \d, b) - /£> ((3, N ± \d, b) (4.1) 



where (0, N ± \d, b) is the excess free energy of the sys- 
tem under periodic boundary conditions when no such 
a diffuse interface exists. Obviously, the helicity mod- 
ulus of the infinite system then simply is T(f3\d,b) = 
liniTv^^oo T (A N±\d, b). Within the isotropic spherical 



model the corresponding result for Y((3\d) is known, sec, 
e.g., [27] 



P?(T\d) = -(K-K c ). 



(4.2) 



The needed information for fcx(P, N±\d, b) is also avail- 
able, see, e.g., [17] 



0fg>(J3,N ± \d, b) 



(d-i) ( b±_ 



(d-l)/2 



x (v d > 2 - v d/2 \ - v d/2 2 V Kd ' 2{n ^ ] } 

x \y P Vb ) v P (27r) d / 2 ^ {n^y/^ y 

(4.3) 

where y p is the solution of the equation 

2 



i rq-d/2) d/2 _, d/2 _, 



E 



K d /2-i{n^/y^) 



(4.4) 



Using Eqs. (3.22) and (4.3), for the finite-size scaling 
behavior of the helicity modulus we obtain 



(3T((3,N ± \d,b) = N ± (d 2) ^ 



j \(d-l)/2 



Xr(x t ), (4.5) 
where the scaling function of the helicity modulus T is 



1 



x r(xt) = ~2< -x t {y-y p ) - - 



y d/2 



lT(-d/2) 



(2TT) d / 2 



2 (4^) d / 2 (?' 
^ [ ' {nVy) d / 2 



<<" 2 -y d P /2 ) 



-v d P /2 J2 



n=i 

Kd/2(ny/%) 



(4.6) 



where y is the solution of Eq. (3.35), y p is the solution 
of Eq. (4.4), and x t is defined in Eq. (3.23). Taking 
into account that when T < T c and N± 3> 1 one has 
y p — > + and y — > + , from Eqs. (4.5) and (4.6) one 
derives, within the spherical model, the behavior of the 
"bulk" helicity modulus in an anisotropic system 



0?(T\d,b) = ±b ± (K-K c ). 



Results for the case d = 3 



(4.7) 



Since d = 3 is of special importance we, similar to what 
we have done for the Casimir force in systems with diffuse 
interface, will present in more details explicit results for 
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the finite-size behavior of the helicity modulus in this 
case. When d = 3 the equations (4.6) and (4.4) simplify 
to 

pT(/3,N ± \d = 3,b) 

= "i'(^)l{r*-»)-ii(* s/2 -^ 2 ) 

V^Li 2 (-e-^ 5 ) - V^Li 2 (e"v^) 



1 



+Li 3 (-e _V5 ) - Li 3 (e"v^) 



and 



x t = 7T VVp + - In [1 - e v^] 



(4.8) 



(4.9) 



respectively. The solution of Eq. (4.9) for periodic 
boundary conditions is 



y/y^ = 2 arcsinh 



1 



(4.10) 



which has to be compared with the corresponding so- 
lution for the antiperiodic boundary conditions, see Eq. 
(3.45). 



from Eq. (4.1) one obtains 



p c T(f3 Cl N ± \d = 3,b) 
2 fJ±_ 



TT 2 N ± V J\\ 



J± 



[Acasimir ^Casimir] 



TT 2 N ± V J, 



S<* (!) 



- ().()<S(i(i4<) V 1 ( j-^j . 



(4.12) 



Taking into account the relation g(T c , L) = 
(m/ft) 2 Y(T, L) between the superfluid density frac- 
tion and the helicity modulus (strictly speaking this is 
valid only for n = 2) one can obtain, within our model, 
an estimation of g{T c , L) at T c . 

The dependence of the scaling function Xy(x t ) is 
shown in Fig. 2. It is easy to show that Xy(x t ) de- 
cays exponentially fast for Xt ~> 1, while for x t <C — 1 
one derives that 



X r (x t ) 



-x t /2. 



(4.13) 



The asymptote of X T for T < T c leads to Eq. (4.7) 
for the behavior of the helicity modulus within the 
anisotropic 0(n) models when n — > oo in the limit 
lim^^oc (3T{(3,N ± \d = 3,b). 




T(T,L) for d = 3. One observes that it is a monotonically _Q^5i"rTTVTV"~7~T~7~7~ i i i i i i i i i i i i i i 

decreasing function of x t . The asymptote of Xy(x t ) for x t <€. "1-5 -1 -0.5 0.5 1 1.5 

— 1 is given in Eq. (4.13). x 



Let us determine the critical value of the finite-size he- 
licity modulus T(/3 C , N±\d = 3, b). Knowing the Casimir 
amplitude for antiperiodic boundary conditions Ac as imir 
(see Eq. (3.50)) and that one under periodic boundary 
conditions [7] (see also [19]) 

A%Ln* = -^C(3) * -0.153051, (4.11) 



FIG. 3: The scaling functions Xcasimir(a;t) of the Casimir 
forces F^ silnil (f3,N ± \d = 3,6) and Fg^tf, N ± \d = 3,6) 
for d — 3. The difference is due to the contributions stemming 
from the helicity modulus. We see that this contribution is 
rather strong and dominates the behavior of the force under 
antiperiodic boundary conditions converting it from attrac- 
tive (under periodic boundary condition) into a repulsive one 
(under antiperiodic boundary conditions). 
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In the current article we studied the behavior of the 
Casimir force and the helicity modulus in anisotropic 
system with a diffuse interface as a function of the tem- 
perature. The interaction along the film is characterized 
via a coupling constant Jm while in the perpendicular to 
the film direction it is J±. We have found that all scal- 
ing functions, including the Casimir amplitudes, depend 
on the ration J±/J\\ and are, thus, nonuniversal. More 
precisely, we have found that the Casimir force in a d- 
dimcnsional system, with 2 < d < 4, can be written in 
the form, sec Eq. (3.38) 



F. 



(a) 



Casimir 



(T,N ± \d,J ± ,J l{ 



= {k B T c ) N± d Xcasimir^tM, J±, J\\), (5.1) 

near the corresponding bulk critical temperature T c , 
where x t is a properly defined temperature dependent 
scaling variable and the nonuniversal scaling function 



-V, „ s!> ,j, (.)', \<l. ./_ ( — 



(d-l)/2 



^Casimir 



(5.2) 

can be related to XcasimiT(xt\d), which is the univer- 
sal scaling function characterizing the corresponding 
isotropic system. The explicit form of Xc a simir(^t|rf), for 
2 < d < 4, is given in Eq. (3.32) and, equivalently, in 
(3.34). Similar relations can be written also for the he- 
licity modulus, see Eq. (4.5) 

r(T,N ± \d,J ± ,J\\) = (k B T c ) N- {d - 2) Xr{x t \d,J ± ,J\\), 

(5.3) 

where, again, the nonuniversal scaling function 
Xr{xt\d,J±,J\\) 



X r (x t \d,J ± ,J n )= (Jpj 



(d-l)/2 



X r (x t \d) (5.4) 



can be related to universal scaling function X?(xt\d) 
characterizing the corresponding isotropic system. The 
explicit form of X?(xt\d), for 2 < d < 4, is given in Eq. 
(4.6). 

From. Eq. (5.2) one obtains, see Eq. (3.39), that 

A C asimir(rf| J±, J\\) = \y~J^ J A CaS imir (ci) ■ (5.5) 

Since, within the spherical model, sec Eq. (3.40), 



(5.6) 



all the relations (5.2), (5.4) and (5.5) are in full confor- 
mity with our general prediction given by Eqs. (1.13) and 
(1.14) which relate quantities of one anisotropic system 
to the corresponding ones in the isotropic system. 

In addition to general expressions pertinent to the case 
2 < d < 4, we have also derived explicit results for the 




case d = 3. The scaling function of the Casimir force 
is given in Eq. (3.41) and, equivalently, in Eq. (3.42). 
The behavior of this function is visualized in Fig. 1. The 
scaling function for the helicity modulus is presented in 
Eq. (4.8) and is depicted in Fig. 2. For the value of the 
Casimir amplitude at d = 3 one has, see Eq. (3.50) [45], 



^Casimir 



(!) 



03) 
6ir 



J± 

Jn 



(5.7) 



while the value of the helicity modulus at T c is, see Eq. 
(4.12), 



PcT(T c ,L) 



Ch 



+ 



7C(3) 



30tt 



■h 



(5.8) 

Let us note that both the Casimir amplitude, as well as 
the Casimir force are positive, i.e. they correspond to a 
repulsion between the plates of the system. Let us stress 
that this effect is solely due to the existence of a diffuse 
interface in the system. We recall that under periodic 
boundary conditions for d = 3 and in the notations of the 
current article the Casimir force under periodic boundary 
conditions is given by the expression [17] 



Li 2 (e"v^) - -Li 3 (e" 




(5.9) 



The comparison between the force under antiperiodic and 
periodic boundary conditions is shown in Fig. 3. We ob- 
serve that the contribution of the helicity energy is so 
strong that the Casimir force converts from being every- 
where attractive (under periodic boundary conditions) 
into everywhere repulsive (under antiperiodic boundary 
conditions). This idea can eventually be used for practi- 
cal purposes when applying some ordering external field 
might cause the border spins, dipoles, etc. to order in 
parallel or in antiparallel way to each other. Of course, 
by changing the degree of helicity the force will pass from 
being attractive through being zero into being repulsive. 
Obviously, it will be interesting to consider such a sce- 
nario in more details by say, studying a system under 
twisted at a given angle boundary conditions. We hope 
to return to this problem in a future work. 
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APPENDIX A: EVALUATION OF U (w, N±\d, b) 

In the current appendix we prove the validity of Eq. 
(3.25) for the behavior of U (w, N ± \d, b) when N ± > 1 
and < w <C 1. Because of the representation (3.6) of 
U (w, N±\d, b) and the asymptotes (3.10) of Sffi (x) one 



divides the region of integration in two subregions - from 
to aNj_ and from aNj_ to infinity, where a is a fixed real 
number such that < a < 1. Let us denote the integral 
over the first region (over " moderate" values of x) by U m 
and let Ui is the integral over the "large" values of x, i.e. 
let 



1 f aNl dxi 
" 2 7 x\ 



exp(— x) — exp < — x 



w — b± ( 1 — cos 



N ± J 



e- x ^I a (xb ± ) + 



nxbj_ 



V xb ± , 



x [e-^iUo^H)] 



d-i 



(Al) 



and 

U t 



1 r 00 ^( e -x 

2 /.at* a; [ 6 



~ xw [e- xb U (xb ]] )] d 1 x 



A ± + 7V ± ^27Vi X 



(A2) 



Obviously £/ = Ui + U m . The evaluation of [// is straight- 
forward. Since x » 1 in calculating U\ one can use the 
large value asymptote of the Bessel function [42] 



h{x) 



exp(x — v jlx) 
\/2wx 



1+ 8~x 



9 - 32u 2 
2\(8x) 2 



with the help of which one directly obtains that 



(A3) 



Ut = -NZ 1 



(d-l)/2 



1 



^. x -(d-l)/2 e -yx 



ab ± /2 x 



( 47r )(d-i)/2 

1 + RW ( 



,(A4) 



where y is defined in Eq. (3.18). Let us deal now with 
the term U m . We divide this term into " bulk-like" contri- 
butions U m .b and "finite-size" contributions U m ^ s where 
U m = U m , h + U m f s with 



^Uxr 
x I 



e -xw [e-'^IoixbA.)] x 
[e-^i/o^ii)]" -1 } (A5) 



and 

C^m,fs 

where 



iN, 



dx 



x yJ2irxbj_ 
x [e- xb U ( xbll )] d -\ 



/ 2N]_ 



«) = u; — 6j_ ( 1 — cos 



(A6) 
(A7) 



It is straightforward to evaluate U m .f s . Due to the repre- 
sentation (3.15), for all x <C N]_ the corresponding con- 
tribution into the integral on the right-hand side of Eq. 
(A6) will be exponentially small. Thus, one again can 
use in (A6) the large-value asymptote (A3) of the Bessel 
function Iq(x) which leads to 



U 



m.fs 



r ab ^ 2 dx 

' Jo x 



b ± \ {d - 1)/2 1 



x d/2 e^ x R^ 



~ ■ ( A8 ) 



It remains now only to deal with the term Z7 mj b- By 
subtracting and adding, up to the linear in x term, the 



asymptote of exp 



xb± 



for small values of 



x one rewrites Eq. (A5) into the form 



U 



m.b 




1+ 2 b± Nf X] ' 



x [e- xb ±I Q (xb ± )] [e- xb U ( x b\\)] 



d-l 



i r aN1 dx 



1 + -b±^2X | i -i 



2 N 



± 



x [e- xb ±I (xb ± )] [e-^Ioixb^ \ (A9) 



It is easy to understand that the integration over small 
values of x in the second line of the above equation will 
provide contributions of the order of 0(AJ 4 ) which we 
will neglect, since we are only interested in contributions 
that are not smaller than 0(N]^ d ), with 2 < d < 4. 
Thus, in this integral one again can use the large value 
asymptote (A3) of the Bessel function Iq(x), which leads 
to 
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y m,b 



aNi 



dx 



[e- xb ±I (xb±)] [e- x6 ii/ (a;6||)] 



d-1 



i 2 raNl 

dxe-™[e-* b ±I (xb ± )]x 



[e-^ii/o^n)]' 



i /&, \ (d ~ 1)/2 i i r ab±/2 dx 



2 x) - e"^ 1 x- d / 2 . 



(A10) 



One can complete the integral in the first and second line of the above equation so that the integration is from to 
oo and to subtract the parts of integration from aN]_ to oo. In the subtracted parts one can again use the large value 
asymptote (A3) of the Bessel function Iq{x). In this way one obtains 



1 7r 2 //) i \ (d ~ 1)/2 11 r°° At 

u^MMb) - -^ w w dHb ) + [ j^ /2 - e- (1 + « 2 x) . 

Jo x 



-d/2 



e-*y l x - d > 2 



(All) 

Expressing from Eq. (3.28) function R(~\x) in terms of R^ + \x) and substituting the so-obtained representation in 
Eq. (A8) one obtains 



(b \ (d ~ 1)/2 ] f 

,fs = ~ NI \h) WW^Jo 



<- ab ^ 2 dx 



x -{d-\)/2 e -yx 



X 



1 + JJ(+) (tt 2 ^ 



1 1 



2 -v/4ttx 



N (d-l)/2 1 f .abx/2 j_ 

" X UJ (47T)( rf - 1 )/ 2 W 



1 + fl<+> U 2 X) - -^L= (1 + n 2 x) 



1 1 f ab ^' 2 dx 



2 \f^K Jo x 



—x- d / 2 e- yx 



(l + tt 2 x) - 



(A12) 



In a similar way, by adding and subtracting the asymptote of 1 + R^ (tt 2 x^ for small values of the arument, one can 
rewrite Ui (see Eq. (A4)) into the form 



(d-l)/2 



i i r b±/2 dx 



2 VAttJq 



(47r )(d-i)/2 | J ab±/2 x 
- d / 2 e-y*(l + n 2 x) ). 



to x -(d-l)/2 e -yx 



1 + R^ U 2 x)---^(l + n 2 x) 



(A13) 



r 



By adding U m . h , J7 m ,f s , and {/ z as given by Eqs. (All), 
(A12) and (A13), respectively, one obtains, after us- 
ing the representation (3.20) for Ud(w\b), as well as the 
fact that Wd{w\b) = dUd{w\b) j du , the final result for 
U (w, N±\d, b) given in Eq. (3.25) in the main text. 



APPENDIX B: DERIVATION OF THE SERIES 
REPRESENTATION OF I(y, d) 

In this appendix we derive the power series represen- 
tation (3.30) of the integral I(y,d) defined in Eq. (3.27). 

First, let us note that using the representation (3.11) 
for the function R^ + '(x), the integral I(y,d) can be de- 



composed as I(y,d) = I^(y,d) + I^(y,d), where 
jW(y,d) = 1 



and 

lM(y,d) 



(47r)( d - 1 )/ 2 

oo /■oo 



Y / dxx-^/ 2 e-y x e-^ n+1 > 



(Bl) 



(4 7r )(d-l)/2 

/•OO 

/ dxx- {d+1 
Jo 



)/2 e -yx 



1 + TT 2 X 

2\/Attx 



(B2) 
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Employing dimensional regularization, the latter integral 
can be done analytically and becomes 



(d-2)/2 



)(47r)d/2 



(B3) 



-4v^T[(l-d)/2]+yr(-d/2)}. 



Introducing the variable y = y — tt 2 , the integral I^(y, d) 
can be written as 



1 



(4 7r )(d-l)/2 

00 /'OO 



^ X -j/x -[it 2 +An 2 n(n+l)]x 



(d+l)/2' 



(B4) 



and upon replacing exp(— yx) by its Taylor series repre- 
sentation the integral ^ becomes 



jW(l/,d) 



OO OO 

i V V b {d) 

(d-l)/2 ™." 



(-y) T 



( 47r )(d-i)/2 Z-j m | 
with the coemcients 



(B5) 



b ( nl(y)= (°° dxx m -^l 2 e-^ 2+ ^ n ^+ 1 ^. (B6) 
Jo 

Again in the sense of dimensional regularization the x- 
integration in the latter equation can be performed to 



give 



(27T) 



Inserting this into Eq. (B5), the n-summation can be 
done analytically leading to 



2 W E 



(-y)"T(m+i^) 



-2m™ J 



-E«S?(-»r 



m=0 



(B8) 



with the coefficients a™ defined in Eq. (3.31). The first 
m-sum in square brackets can also be done analytically 
and we obtain 

jNfo.d) = (4 7 r)( 1 - d )/ 2 r[(l-rf)/2](y + 7 r 2 )( d - 1 )/ 2 

OO 

-^- 1)/2 E^(-^) m (B9) 

m=0 

If we now add up /M (y, d) and 1^ (y, d) we arrive at the 
power series representation (3.30) of I(y,d) given in the 
main text. Note that no terms being nonanalytic with 
respect to y are present, and furthermore that the radius 
of convergence of the expansion is \y\ = \y — n 2 \ < A-k 2 . 
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